Evolution of speckle during spinodal decomposition 
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Time-dependent properties of the speckled intensity patterns created by scattering coherent radi- 
ation from materials undergoing spinodal decomposition are investigated by numerical integration 
of the Cahn-Hilliard-Cook equation. For binary systems which obey a local conservation law, the 
characteristic domain size is known to grow in time r as -R = [-Br]" with n = 1/3, where B is 
a constant. The intensities of individual speckles are found to be nonstationary, persistent time 
■ series. The two-time intensity covariance at wave vector k can be collapsed onto a scaling function 

O ' Cav{5t,t), where &t — k^^"B\T2 — nj and t = k^^"B{Ti + T2)/2. Both analytically and numerically, 

the covariance is found to depend on St only through St/t in the small-? limit and St/t^~" in the 
large-f limit, consistent with a simple theory of moving interfaces that applies to any universality 
class described by a scalar order parameter. The speckle-intensity covariance is numerically demon- 
ic , strated to be equal to the square of the two-time structure factor of the scattering material, for 
^ ' which an analytic scaling function is obtained for large t. In addition, the two-time, two-point order- 
parameter correlation function is found to scale as C (^r/ ^B" y^rpM-rf^^ , ti/t2^ , even for quite 

a large distances r. The asymptotic power-law exponent for the autocorrelation function is found to 
be A « 4.47, violating an upper bound conjectured by Fisher and Huse. 



PACS; 64.75.-(-g, 64.60.Cn, 61.10.Dp, 05.70.Ln 



I. INTRODUCTION 

(N 

> 

' One of the most common ways to obtain information about the spatial structure of a material involves the scattering 
I of radiation, such as electrons, neutrons, or photons. The scattering is caused by inhomogeneities in the scattering 
CO . cross section of the material. Depending on the nature of the incident radiation, this scattering cross section is 
1^ ' given by the local value of some property of the material, such as its density or magnetic polarization. In the first 
Born approximation the amplitude of the scattered wave is given by the Fourier transform of that property, and the 
wave vector k is proportional to the momentum transfer as a coherent wave of the probe radiation is scattered. For 
conventional radiation sources, the wave is only coherent over a small volume and different parts of the sample scatter 
independently of each other. As a result, only the ensemble average of the scattering intensity can be observed. 
On the other hand, when the coherence volume of the wave is sufhciently large to encompass the whole scattering 
I volume, the unaveraged Fourier transform of the specific geometric structure of the material is observed. If the 
'"O inhomogeneous scattering property of a particular sample at position r and time t is represented as the scalar field 
^ '(/'(r, r), and its Fourier transform at wave vector k as i/'(k, r), then the observed scattering intensity is 

O 



/(k,r)= ij{k,r) ' , (1) 



where proportionality constants have been ignored. The domain pattern of the order parameter ^'(r, t) from one of 
^ I our simulations of spinodal decomposition in a two-phase material is shown in Fig. |l| with the shading differences 
■ - ' corresponding to differences in the local scattering cross sections. Figure |^ shows the central section of the scattering 
intensity corresponding to Fig. |^, with darker shades indicating higher intensities. The speckled appearance of the 
scattering intensity is a characteristic feature in the scattering of coherent radiation. 

Since the speckle pattern depends on the specific configuration of the material, it changes as the domain structure 
evolves over time. The result is that the individual speckles fluctuate around their time-dependent averages. The 
intensity fluctuations are almost uncorrelated in k-space, but the evolution of the speckle pattern gives rise to nontrivial 
two-time correlations in the scattering intensity at individual wave vectors. 
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The fluctuation of a single speckle around the average intensity is shown in Fig. ||. It has been normalized at each 
time by the corresponding average intensity. Also shown in Fig. |^ are Brownian noise fluctuations, whose two-time 
covariance is exponential. That noise is constructed to have the same single-time probability density and characteristic 
time scale as the speckle. There is a clear qualitative difference in time correlations between the two time series. The 
changes in the Brownian function are relatively large and the fluctuations shorter lived than those of the speckle 
intensities produced by the phase-ordering simulation. This highlights a property of the speckle intensities called 
persistence [0. While the Chan ges in a Brownian time series are independent, the changes in a persistent time series 
are correlated with each other. 

In the past, photon correlation experiments using laser light at wavelengths ranging from the infrared to the 
ultraviolet have used the two-time correlations of individual speckles to study fluctuations in many materials . 
New high-brilliance synchrotron photon sources have made coherent x-ray experiments feasible. The comparative 
advantages of x rays over lasers for many systems include penetration of optically opaque materials and the ability to 
investigate shorter length scales. Coherent x-ray intensity fluctuation spectroscopy (XIFS) experiments have already 
been performed on a number of materials . Recently, the two-time correlations in speckle intensity were measured 
by XIFS for a sodium borosilicate glass undergoing spinodal decomposition 

Our intent with this paper is to make a theoretical study of how speckle experiments can be exploited to investigate 
the relaxation of homogeneous materials towards a new, heterogeneous equilibrium after a rapid change in a thermo- 
dynamic parameter such as temperature. Previously P,p^ this was done for materials subject to phase ordering, in 
which case ip{r, t) is not constrained by any conservation laws (often called model A Here we analyze scattering 

from materials subject to a local conservation constraint (model B such as the process of phase separation in 

a model binary mixture or alloy. The dynamics of the phase separation that occurs after a high-temperature ho- 
mogeneous mixture is suddenly brought into the low-temperature regime of two-phase coexistence is a fundamental 
statistical-mechanical process that has important consequences for multi-phase materials. Two examples are metallic 
alloys jTT| and Vycor glass p^ . 

The scaling properties of correlations in the model system and their relationship to speckle correlations are discussed 
in Sec. Although our focus is on model B, many of our results are sufficiently general that the also apply to other 
universality classes characterized by a scalar order parameter, such as model A. Details of the numerical simulations 
are given in Sec. IB. In Sec. IV the numerical results are presented and compared with the theoretical scaling relations 
from Sec. Our conclusions are summarized in Sec. 



II. SCALING THEORY 



When the volume fractions of the two components are equal in a binary mixture, the homogeneous high-temperature 
phase is unstable to long-wavelength fluctuations after the quench. As a consequence, an interpenetrating pattern 
of domains rich in one or the other of the two components develops immediately after the quench. (See Fig. 1.) 
This intertwined structure distinguishes spinodal decomposition from many other phase-separation processes. At 
later times, the domains grow as the system evolves towards an equilibrium state with the smallest possible total 
interfacial area compatible with local conservation of the volume fractions of the constituent phases. The patterns at 
different times are statistically similar, except for the characteristic length scale of the domains. The description of 
the late-time regime of phase separation is greatly simplified by this dynamic scaling behavior. 



A. Scaling of one-time functions 

Usually, the average size and shape of the growing domains can be measured using the structure factor, which is 
proportional to the average intensity measured in scattering experiments. For isotropic materials the structure factor 
depends only on the magnitude of the wave vector, 

5(fc,T) = (/(k,r)) , (2) 

where the angular brackets represent the ensemble average. For locally conserved dynamics, the structure factor 
vanishes at both high and low wave vectors. The intervening maximum occurs at fcmax(''"), which is proportional to 
the inverse of the characteristic domain size. Since the domains grow, the ring of brightest scattering contracts as the 
system evolves toward equilibrium. 

In the scaling regime, the growth of this characteristic length obeys a power law, R{t) — [i?r]", where B \s a. 
constant. Such power-law growth is found in many different phase-ordering processes, and the exponent is usually 
independent of the specific properties of the material. Instead it depends on general features of the phase ordering. 
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such as the presence or absence of local conservation. This dynamical universality picture has been found to be widely 
applicable, and a group of diverse processes that show common scaling behavior are termed a universality class 
For situations described by a scalar order parameter, n = 1/2 when the order parameter is not conserved and n = 1/3 
when it is locally conserved. 

Dynamical scaling and dynamical universality can be combined to give a single, time-independent description 
of the statistical structure of phase-separating materials. Since R{t) — [St]", a dimensionless scaled wave vector 
q = kR{T) = fc[Sr]" can be defined and the structure factor can be collapsed onto a master curve 

F[{q) = [BTr^''S{k,T) , (3) 

where d is the spatial dimension of the system. Equivalently, a scaled time t — q^/"^ — [kR{T)]^^'" — k^^^Br gives the 
master curve 

F,{t)^k'S{k,T) , (4) 

where Fi{t) = q'^F[{t'^). Incoherent scattering experiments are frequently used to measure the average scattering 
associated with the structure factor. 

B. Scaling of two-time functions 

For individual speckles, time-dependent fluctuations around the average intensity are characterized by the two-time 
intensity covariance function. 



Covk(k,ri,r2) - ^/(k,ri)/(k,T2)^ - ^/(k, n)^ ^/(k, T2) 
J^(k,Ti)^*(k,ri)V'(k,r2)V'*(k,T2)) 
-(^(k,ri)V;*(k,ri))(v;(k,r2)^*(k,T2)) . (5) 

Following the treatment in Rcf. [p^ , we anticipate that the joint probability density of the complex variables '!/'(k, ti) 
and '(/'(k, T2) is Gaussian. The Gaussian nature of the joint probability density, and thus of the marginal probability 
density for each '0(k, r), is essentially a consequence of the central-limit theorem. It may or may not hold, regardless 
of the probability density of the order parameter in real space V'(r,T). Specifically, it does not depend on whether 
'0(r, r) can be derived from an underlying Gaussian auxiliary field, an assumption which is customarily referred to as 
"the Gaussian approximation." In this paper we are only concerned with the former approximation. 

It is possible to construct a set a variables whose joint probability density is not Gaussian even though the marginal 
probability densities are |jl^, so the additional constraint that the joint probability density be Gaussian is an important 
one. The essential aspect of this assumption is that the fourth moment which occurs in Eq. reduces to products 
of second moments according to Wick's theorem. For this reason we refer to the approximation as the "Gaussian 
decoupling approximation." For k 7^ the covariance becomes equal to the square of the two-time structure factor 
0, 

Covk(k,ri,r2) = 52(k,ri,T2). (6) 
The two-time structure factor is defined as 

5(k,ri,r2) = (^(k,ri)V;*(k,r2)) , (7) 
but it is also proportional to the Fourier transform of the two-point, two-time order-parameter correlation function 

5(k,Ti,r2) = j rfre* '-C(r,Ti,r2) , (8) 

where the correlation function is defined by 

C(r,Ti,r2) = (V(0,ri)V(r,T2)) . (9) 

Similar to the case for the one-time structure factor, dynamical scaling also applies to the two-time structure factor. 
For isotropic media, the latter can be related to a master function of the two scaled times. 
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F2{h,t2) = k''S{k,Ti,T2) 



(10) 



Since it involves the Fourier transform of the material at two different times, the two-time structure factor cannot be 
measured directly in a scattering experiment. However, its relationship to the two-time intensity covariance indicates 
that it can be inferred from the fluctuations of speckles around their average intensities. In analogy with the scaling 
of the two-time structure factor, a scaled speckle covariance 

C0vih,t2) = fc2'^C0Vk(fc,Ti,T2) (11) 

can be defined. Then 

Cov{h,t2) = Fi{ti,t2) , (12) 

provided the relationship between the speckle covariance and the two-time structure factor, Eq. (|^), is valid. 

The asymptotic properties of the two-time structure factor can be found using quite general arguments that should 
apply for many universality classes. They are most concretely expressed by focusing on the geometry of the moving 
interfaces, a schematic view of which is provided in Fig. ^. The top and bottom of the figure remain inside the domains 
of the two different phases, respectively, while the shaded region changes phase between ti and T2. The edges of the 
shaded region represent the interface between the domains at ti and T2; the interface at T2 is flatter. Roughly halfway 
between these interfaces, the dashed line represents the "mean interface," which defines a mean domain size, [-Bf]", 
where f — {t2 + ti)/2. Lengths are scaled by this mean domain size. Using this specific choice of scaling length, 
along with the ansatz that the times appear as the ratio T2 /ti , the relationship between the two-time order-parameter 
correlation fimction and its master curve can be written 

C{r,n,T2)=H{r/[Bfr,ST/f) , (13) 

where St — \t2 — ti\. 

The master curve for the two-time structure factor can be found using Eq. (H). When the correlation function 
is isotropic, this can be reduced to a radial integral that depends on the scaled-time variables St ~ \t2 — ti\ and 
i= (ti +<2)/2; 

F2((5t,t) = (2^)'*/2t"W2+i) / duu'^^^Jd/2-i{ut")H{u,St/t) , (14) 

where Ji, is a Bessel function of the first kind of order v. 

The small-i and large-? behaviors of F2 can both be found using elements of the picture described here. 



1. Two-time scaling for small t 
The small-i behavior of Eq. (nj) is isolated by Taylor expanding the Bessel function. 



F2 {St, t) = 27r''/2f V(-l)^"P'y (St/t) , (15) 



J=0 

where h2j{St/t) is defined through the integration over the scaled distance 

duu^^+'^-^H(u,St/i) 

'■''W- ° 4^,!r(, + d/2) - <") 

At small wave vector the structure factor can be shown to obey a power law S oc k^^^\ or equivalently F2 oc f "■(/'(")+'^). 
The value is /3(l/2) = for model A |jl^ and believed to be /3(l/3) = 4 for model B For model B, this indicates 
that ho{0) = ^.2(0) — 0, and, consequently, that the leading term is the one containing hn{St/t). The Taylor expansion 
is well approximated by its leading term when that term is much larger than the one containing ft,g. This condition 
is equivalent to 

-2„ fd \ C duu'^+^H{u,St/t) , ^ 

where the ratio of the integrals is a function of Stjt. For t small enough, then, the scaling of the two-time structure 
factor depends on St only through Stjt. The analogous result for model A (which follows with ft.o(0) > 0) was already 
obtained explicitly in Ref. [|l^ for the Yeung-Jasnow (YJ) ||l^ correlation function. This small-t behavior agrees with 
our numerical results for both model A |l^ and model B ^ee Sec. IV). 
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2. Two-time scaling for large t 

When Eq. Q is not satisfied, a large number of terms is required in Eq. ( |l5|) to obtain an accurate estimate of 
F2{St,t). Since these terms contain products of powers of i and St/t, F2{St,t) should not be expected to depend on 
St only through 5t/t for large t. In Ref. this was demonstrated explicitly for model A using the YJ analytic 
approximation for the correlation function. The result presented here is based on simple geometric arguments, which 
should be valid for all universality classes with a scalar order parameter, including both model A and model B. 

At large wave vectors the one-time structure factor for a two-phase system with sharp, randomly oriented interfaces 
is well described by Porod's law jl^, S oc Here the proportionality constant includes the interface area 

per unit volume (the specific surface) at time r, A{t) oc \/R(t). In this limit the scattering probes the correlation 
function at lengths r <C R{t), where the normalized equal-time correlation function decreases linearly with scaled 
distance ^(r)r cx r/i?(r) 

Here Ci is a nonuniversal constant. The mean-square order parameter {iP'^{t)) is independent of r if the two equilib- 
rium values of ip are equal in magnitude or the volume fractions of the two phases are constant, as they are in Model 
B, and we will simply denote it ("0^). 

The behavior for small, but nonzero, time differences can be deduced from the moving-interface model. The 
normalized two-time autocorrelation function C{0,ti,T2)/ {'ip'^) is reduced from unity by an amount proportional to 
4>A, the fraction of the system volume which is occupied by different phases at ti and T2, i.e. the shaded volume in 
Fig. ^. In the sharp-interface picture used here, this is an exact result, obtained by simple probability arguments. 
The proportionality constant is easily calculated but not very useful. 

For small time differences in a moving-interface model, one can construct a coordinate system based on the mean 
interfacial position, represented by the dashed curve in Fig. |^, and a coordinate locally perpendicular to it, A. The 
specific surface of the mean interface is approximately A w -I- A{t2)]/2, and the average distance between the 

other interfaces is (|A|). In this approximation (j>/^ « A (|A|). Using the scaling assumption that R{t) ~ [Bt]" cx 
1/A{t) is the only length scale characterizing the system at time r, we obtain 



(|A|)(x|i?(r2)-i?(n)|- '^^^^^ 



dr 



5t = nB''J^=nB''C . (19) 



The dimensionless proportionality constant between (|A|) and \R{t2) — R{ti)\ is expected to depend on the particular 
dynamic model and the dimensionality d. The relation \R{t2) — R{ti)\ « nB"(^ is a Taylor expansion valid for small 
(. These results can be combined to give 

g(0,n,T2) _ _ nC,C2C 

(02) ^« ' ^^^) 

where C2 contains the proportionality constants from Eq. (p^. 

On length scales (|A|) ^ r ^ 1/A the interface should appear sharp, so the asymptotic correlation function should 
be described using the mean domain size, i.e. by Eq. ( p^ ) with r~" replaced by f ~". The crossover between the 
small-r and intermediate-r limits can be described in terms of a scaling function G{x) as 



(02) V"C'2S"C 

with the asymptotic behavior 

Thus, C(r, (5t, f) depends on r only through the dimensionless scaling combination r/[nC2-B"C]- As seen from the 
derivation, this result is general and applies to any moving-interface model. 

The two-time structure factor in the large-fc limit is obtained from Eq. ( pT| ) using the same formula for the d- 
dimensional Fourier transform of an isotropic function as in Eq. (p^, 

5(fc, n, r,) = -fc-'^(02) (2^)^/2 duu'l^ G (^^^^) , (22) 
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where we have removed from C(r, St, f) the constant term which only contributes to a 5 function at = 0. Defining 
the new dimensionless scaling variable z = kB'^( = St/t^~'^, we obtain the asymptotic scaling form for the two-time 
structure factor in the large-? limit: 



nC2Z 



(23) 



This asymptotic scaling function depends on k, ti, and T2 only through z. If the correlation function vanishes 
sufficiently smoothly for large r, standard expansions of Fourier transforms pC|j2l[ | indicate that the asymptotic large- 
k behavior of the structure factor is determined completely by the real-space behavior of the correlation function for 
small r. For z = it reduces to the single-time scaling function ||20t] 



$d(0) 



t"f2(0,^") 



2<i7^(rf-l)/2p 



1 



(24) 



which is consistent with the amplitude of the large-fc Porod tail of the structure factor ($2(0) — 2tt and $3(0) — Stt). 
As z increases, the integral in Eq. (^3|) tends towards the Fourier transform of a constant, which vanishes for all 
nonzero k. Thus, ^d{z) decreases from its maximum value at z = towards zero as z increases. 

In Ref. the scaling behavior expressed in Eq. (|2^) was explicitly obtained for model A, but a nonrigorous 
scaling argument suggested it should hold for other dynamical models as well. This was recently confirmed by XIFS 
experiments in a sodium borosilicate glass undergoing spinodal decomposition [|| . This is a two-phase system with a 
locally conserved scalar order parameter appropriately described by model B. 

To predict the form of ^ d(z) in further detail it is necessary to know how G{x) interpolates between the limiting 
behaviors given by Eq. (21b). A detailed calculation based on an extension of an approach used by Tomita to calculate 



equal-time correlation functions |19| , ^0| , indicates that G{x) — x (x x ^ as x —> 00 and G{x) — 1 oc a:^ as cc ^ for 
moving-interface models 
and the substitution 4 



A convenient analytic form with this behavior is G{x) 
= (wy)^, the integral I{v) = ^d{v/nG2) is converted to 



Choosing v = nG2Z 



I{V) = -{2^Y'^ / du u'^'^Jd/2-l{u) \Jv? + W2 



(2^)./2^./2+2 ^ ~ y2 ^^/^y^'-' J,/,_^ 



(25) 



Asymptotic convergence is guaranteed by the discussion following Eq. (E4 
evaluating a generating function for the desired integral as an Abel limit 



and it is possible to use the approach of 
123|,P| 



I{v) = -(2^)'^/2^;'^/2+2 Hni — dy e""^ ( Vy^) ^<i/2-i [vVy^) • 
This integral can be found in tables |^^, and differentiation followed by taking the limit leads to 

(nC2z)('^+i)/'i^(d+i)/2 {nG2z) 



'i>d{z) = MO)- 



2(<i-i)/2r((d+ i)/2) 



(26) 



(27) 



where Ki, is a modified Bessel function of the second kind. For d — 3 this can be numerically evaluated without 
difficulty, while for d = 2 the expression simplifies further to |26| 



$2(2) = 2tt {1 + nG2z) e- 



(28) 



The same result was obtained in Ref. iQ, using the full YJ ||l^ correlation function for model A (see Appendix). 
The advantage of the present approach is that it demonstrates that ^d{z) is manifestly independent of the large-r 
behavior of the correlation function, requiring only that it converges sufficiently smoothly to ensure the equality of 
the Abel limit and the integral. 



III. NUMERICAL PROCEDURE 



In phase-separating systems no material is created or destroyed. The order parameter is locally conserved, and 
evolution occurs by diffusion along chemical-potential gradients. The Cahn-Hilliard-Cook model [P^,p8l is a convenient 
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description for the dynamics of a conserved scalar order parameter. The thermodynamics of the system are described 
by the Ginzburg-Landau- Wilson free energy ||Tl[] , 



T[i^{v,T)] = / dv 



^V'(r,T) + Jv^4(r,r) + ^|VV'(r,T)p 



(29) 



For a > 0, the first two terms of the integrand create a bistable local potential energy. The last term represents the 
surface tension between domains in which has opposite sign. The dynamics are implemented using the Langevin 
equation 

^.MV^^^t^+,(r,r). (30) 

The first term on the right-hand side represents deterministic relaxation of the chemical potential; the Laplacian 
expresses the local conservation constraint. Processes operating at shorter time and length scales are considered 
thermal noise and are modeled by the random variable ry, whose intensity is given by a fluctuation-dissipation theorem 
^t| . We neglect rj because the most important sources of noise here are the initial conditions, which give the random 
domain morphology at early times. 

For the symmetric mixtures considered here, the parameters can be eliminated by appropriately normalizing the 
time, length, and concentration scales ||l^,|2^,^. Using the same names for the new variables, the resulting dynamical 
equation is 

^^ = -iv^[(l + V^)V'(r,r)-V'3(r,r)] . (31) 

In particular, this yields ip — ±1 for the equilibrium values of the order parameter. 

We have simulated this model on a square lattice, with Ar=l and — Ly = L = 512, using a simple Euler 
integration scheme with At = 0.05. The initial condition was implemented by choosing random values uniformly 
distributed between ±0.1 for each point on the lattice. Measurements were made every 50 time units out to a maximum 
of r = 4000. As can be seen from Fig. ^ the domains are much smaller than the system size, even at this latest time. 
Finite-size effects are therefore not expected to affect the order-parameter dynamics. Indeed, no deviations from the 
expected behavior of the characteristic length are observed in the simulations. 

The other numerical procedures are identical to those described in Ref. JlOt . The Laplacian was implemented using 
eight-neighbor discretization ||3l| , |3^ , and the magnitude of the wave vector fc(k) was defined in a manner consistent 
with that Laplacian. In addition, the Fourier transform of the hardened order parameter, sgn[7/;(r, r)], was used to 
minimize the effect of finite interface width on the scattering intensity. 



IV. RESULTS 

The Gaussian decoupling approximation that leads to Eq. ^ has not been directly justified for phase-ordering 
systems. Gaussian variables are often assumed in the context of the central limit theorem, where it can be argued 
that many independent random variables with finite variance contribute additively. Since R{t) represents the average 
domain size, at any given time there are on the order of {L/R{t)Y independent domains in a system of edge-length L. 
This number can be quite large, but since the domains interact as they grow, with material diffusing along chemical- 
potential gradients, the applicability of the central limit theorem is not obvious. In our previous study of a system 
with nonconserved order parameter | [lO| ], the decoupling was justified a posteriori from the numerical results. Similar 
justification is presented here for numerical integration of the equation of motion for the conserved order parameter. 
When discussing the numerical results, we have taken B = 1 for convenience. 

One way to test the validity of the assumption at the single-time level is looking at the probability density of 
the speckle intensities. Namely, Eq. (^) with ri = T2 is satisfied if the normalized scattering intensity, s(k, t) — 
/(k, T)/S'(k, t), has an exponential probability density, 

P(,s)=exp(-,s), (32) 

that is independent of (k, r) . A single histogram was constructed for the normalized intensity for all times and for 
0.08 < A;(k) < 0.75 |^^. The intensity at a particular (fc, r) was normalized by S{k, t) estimated by circular averaging 
only for the same quench experiment. The lower k cutoff was chosen so that at least 20 speckles contributed to 
that circular average, while the upper cutoff was set using results for ti ^ T2 (discussed below). Nearly 10* samples 



7 



contribute to the histogram in Fig. ||, and the agreement with the exponential form is remarkable. Other histograms 
were constructed to check for k and r dependence, but are not shown here. Wave vector dependence was investigated 
by sampling from narrow rings of constant k at all times, and time dependence was checked by sampling all k 
between the cutoffs at single times. The probability density for the normalized intensity does not appear to depend 
on these variables. A similar check of the Gaussian decoupling assumption was made by Shinozaki and Oono | |3^ for 
a cell-dynamical simulation with conserved order parameter. 

The equality between the two-time structure factor squared and the two-time intensity covariance implied by 
the Gaussian decoupling assumption can also be tested directly. These two quantities can be estimated from the 
simulations using Eq. (|^) and Eq. respectively, with circular averaging over k in addition to the averaging over 
the quenches. Since the real-space correlation function is an even function of r, S{\s.,ti,T2) is real valued. We have 
verified that the imaginary part is zero to within the accuracy of our results. The equality is tested graphically in 
Fig. H as a function of ti = k^Ti for (a) ti = 100, T2 = 200 and (b) ti = 2000, T2 = 4000. The intensity covariance 
and the squared structure factor are equal within our numerical accuracy for sufficiently small k. However, for z ^ 4.5 
(marked by arrows in both parts of Fig. ||) the covariance becomes larger than the squared structure factor, indicating 
the gradual breakdown of the Gaussian decoupling approximation for large z in these simulations. A systematic 
investigation of the range of validity of the Gaussian decoupling approximation, including the possible system-size 
dependence of the number of independent contributions to the Fourier transform, is left for future study. In this paper 
all quantitative scaling results are based only on data for A:(k) < 0.75 |^^, where the squared two-time structure factor 
appears to equal the speckle covariance for all times considered. 

The scaling function ^2(^1, t2), which describes both the two-time structure factor and the related speckle covariance, 
depends on two rescaled times. A normalized analog of the scaled covariance, the correlation function, can be defined 
as 



7(fc,ri)/(fc,r2) 

Corr(ti,t2) = ^. ^-1- (33) 

/(fc,ri))(/(fc,T2) 



Since the Gaussian decoupling approximation for scattering has been verified for ri = T2, this quantity is unity by 
construction for ti = t2- Thus the contour plot of Corr(ti,t2), Fig. ^, is more illustrative than one for Cov(ii,t2). 
The correlation function decays as one moves away from the ti — t2 diagonal. Given the symmetry under exchange 
of ti and t2, t and 5t are more natural variables for the scaling functions Corr and Gov. In particular, t measures 
distance (in units of scaled time) along the diagonal, while 5t measures the distance perpendicular to the diagonal. In 
Fig. it is apparent that the speckle intensity stays correlated for larger 5t as tis increased. As noted in Ref. all 
the quantities used here are readily obtained in experiments, and the method of data analysis should be well suited 
for experimental analysis as well. 

The asymptotic scaling predicted for large and small i can also be tested. First, the characteristic time difference 
5U. defined by Cov((5tc,f) = iGov(0,i) can be found for fixed ranges of t. The measured 5tc^ found by linear 
interpolation, is presented as a function of t on a log-log scale in Fig. |8[ The results show obvious power law behavior 
at both small and large t. A least-squares fit for i < 50 gives an exponent of 1.02 ± 0.02, which agrees with our 



expectation that in this regime the two-time structure factor scales with 6t/t (Sec. II Bl). Least-squares fitting for 
t > 300 gives an expone nt of .66 ±0.01 in excellent agreement with the theoretical prediction oil — n — 2/3 expected 
from z = 5t/t^~^'^ (Sec. pi B 2 ). The fit also indicates that the characteristic t for crossover between these regimes is 



ic ~ 60, which corresponds to the shoulder in the structure factor (see Figs. 6 and 10). 

The scaling function ^d{z) associated with the speckle covariance at large t is shown in Fig. |9|(a) for several t > 700. 
Here the normalization was determined from the asymptotic data where it appears that (2t:CiY ~ 32. The collapse 
of the data is quite good, and the agreement with the theoretical scaling form $2(2) from Eq. ( p8|) is also excellent for 
z < 4.5. These results are shown on a linear-log scale in Fig. ||(b), along with the scaled two-time structure factor data 
for the same i (drawn as filled symbols). The fitting parameter nC2 ~ 0.62 was determined from a nonlinear fit to the 
two-time structure factor data. Fitting to the covariance data did not significantly increase the agreement with $2(^)- 
This theoretical scaling form gives much better overall agreement with the data than more conventional choices, such 
as a Gaussian or Lorentzian. Beyond z « 4 the intensity covariance is quite small. However, as z increases further it 
becomes larger than the squared structure factor, which remains reasonably well described by ^2(2) for all z studied. 
This deviation signals the breakdown of the Gaussian decoupling approximation, and is the same deviation seen at 
large ti in Fig. ^. 

Measurements that include intensities only for two measurement times, Ti and T2, correspond to rays of slope ti/t2 
in the (ii,i2) plane. The wave vector then serves as the parameter indicating position along the ray, and the values 
of the intensity covariance along the ti = t2 ray are the structure factor squared. Figure |lO|(a) shows the simulation 
estimates for the scaled structure factor Fi (t) for several different times on a log- log scale. The collapse for different 
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T is quite reasonable. The two straight hnes are provided as a reference to the expected asymptotic behavior. The 
dashed Une in Fig. |l^(a) indicat es the slope corresponding to the power law S cx /c''("\ or Fi cx ^jth /? = 4 



for model B as discussed in Sec. [IB 1. The solid line corresponds to the Porod's law result S oc fc^'^''+^^ (or Fi cx 
at high wave vectors, which is an important part of the moving-interface model associated with the asymptotic scaling 
form ^d{z)- The simulated structure factor reproduces both behaviors. The inset is the representation of the master 
curve F{ in terms of the scaled wave vector q = fc[i?T]^/^. The deviation from the expected power law at small q has 
been observed in other simulations and is expected to vanish as r — > 00 JlBf . 

In the case of nonconserved order parameter, the analytic theory of Ohta, Jasnow, and Kawasaki (OJK) [p^ , 35| 
gives excellent quantitative agreement for the one-time structure factor [^,0 . The YJ extension of that theory 1 1^ 
to two-time correlations gives good agreement for the intensity covariance except at very large St jl^]. No comparably 
successful theory exists for the conserved order-parameter case considered here. Approximate forms for the one-time 
structure factor (none of which has the expected fc^ behavior at small k) have been proposed by Ohta and Nozaki , 
and by Tomita |^^, Yeung, Oono, and Shinozaki |Q. Of these, we find the best agreement for the latter result. The 
heavy solid curve in Fig. |l^(a) represents this approximate structure factor as obtained from the differential-equation 
formulation in fc-space found in Ref. [^o|| . This theory extends a model of interface motion to phase separation and 
recovers n = 1/3. Here it has been fit numerically to the simulation data using one adjustable parameter. While 
it gives reasonable agreement at small distances, it fails at larger distances as seen by the deviations at small t in 
Fig. |l^(a). There is also noticeable disagreement around the shoulder at high t. A better analytic model for the 
locally-conserved dynamics of phase separation is clearly still needed. 

The scaling of the two-time structure factor along the ray ti/t2 = 1/2 is shown in Fig. |l^(b) for several ri. The 
data collapse is again quite good. The largest change from the one-time structure factor shown in Fig. [To|(a) is the 
disappearance of the Porod tail at large wave vector, which is consistent with the idea that moving interfaces are 
controlling the two-time properties of the system. 

It is also interesting to consider scaling for the two-time order-paramete r corre lation function for large r, where it is 



not linear in r and the scattering is not governed by Porod's law. In Sec. IIB 1, we argued that the scaling variables 
should be r/[Bf]" and St/t, and this gives quite good agreement with the simulations. However, an equivalent set 
of scaled variables gives what we believe is a more natural form for the correlation function. It is inspired by the 
YJ correlation function appropriate to model A Eq. (^6|), where scaled lengths enter as r^/[i?^(Ti) -I- R'^{t2)]. 

(This is a consequence of the assumption of a Gaussian auxiliary field.) The simulated two-time order-parameter 
correlation function for model B is well scaled out to quite large r using 

C(r,ri,T2) = C ( , ,-] . (34) 

The results for fixed ti and several T2 are presented in Fig. pT|(a). The zeros of the correlation function are approx- 

/ 2/3 2/3 

imately stationary with respect to a; = r 

/V^i +^2 for the range of times shown here. The scaling becomes less 
good for much smaller values of ti/t2. As the ratio ti/t2 decreases the amplitude of the oscillations also decreases. 
The collapse of C(r, ri,T2) onto C{x,Ti/t2) for several pairs of ti and T2 is shown in Fig. l^(b). The simulation 
results appear to converge for later times, and the collapse is quite good for ri > 500. The solid curve represents the 
simulation estimate of the master curve for this ratio of times with T2 — 4000. 

One last property of phase-separating materials to consider is the autocorrelation of the order parameter, C(0, ti, r2), 
for T2 ^ Ti. In this regime (the opposite of the small time-difference regime emphasized elsewhere in this paper) 
Fisher and Huse |^ have argued that the autocorrelation should be described by a power-law exponent A as 

C(0,Ti,r2)(x (^||^)' ■ 

The YJ correlation function for model A (see Appendix) predicts A = d/2. This is a weakness of the YJ theory; the 
value observed in experiments and simulations on model A indicate A is clearly larger |l0[ . For model B, Yeung, Rao, 
and Desai (YRD) |4^] found A w 4 for ti in the scaling regime from a numerical simulation similar to the one presented 
here. The autocorrelation for the present simulations is shown in Fig. where the results have been normalized 
by C(0,ri,ri) to show the collapse for different ti. The inset shows the values of A estimated from the derivative of 
C(0,Ti,r2) as a function of T2/T1. For the two earliest ti, A appears to have converged to its asymptotic value. A 
least-squares fit to the latest times {t2 > 2000) for ri = 50 gives A = 4.47 ± 0.03. Both our result and that of YRD 
violate the upper bound A < d conjectured by Fisher and Huse for phase-ordering systems They are consistent 
with the lower bound A > {d/2) + 2 found for spinodal decomposition in d > 2 by YRD |4|]. Lee and Rutenberg ]4^ ] 
found the same bound for small, but nonvanishing, minority-phase volumes. They point out that correlations in the 
minority domains accelerate the decay of the autocorrelation. 
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V. CONCLUSIONS 



The scattering intensity for two-dimensional systems undergoing spinodal decomposition has been investigated 
using numerical integration of the Cahn-Hilliard-Cook equation. The results indicate that the two-time structure 
factor can be measured experimentally via the speckle-intensity covariance at k ^ 0. The connection between the two 
quantities is established directly by comparing numerical estimates of both and indirectly by verifying the exponential 
probability distribution of the scattering intensity. For the present simulations, the equality between the two-time 
structure factor and the speckle-intensity covariance is observed to break down at z ^ 4.5. These properties were also 
observed in previous simulations for non-conserved dynamics , and they should be common to any heterogeneous 
material with a large number of independent domains. 

An important conclusion of our study is that the scaling behavior of the speckle-intensity covariance at both large 
and small t can be obtained from geometric arguments based on a simple moving-interface picture. These scaling 
results depend on the dynamics only implicitly through the values of the dynamic exponent n and the exponent 
which gives the asymptotic small-fc behavior of the one-time structure factor. Specifically, we have argued that the 
asymptotic scaling of the two-time structure factor for large t is in terms of z — 6t/t^~^, with 1 — n — 2/3 for 
spinodal decomposition. For small i, on the other hand, the appropriate scaling variable is 6t/t. These behaviors 
were recently confirmed by an XIFS experiment Q of spinodal decomposition and demonstrate that experimentally 
observable aspects of both the short and long length- and time-scale regimes of phase separation can be described by 
relatively simple models of interface motion. A simple form motivated by the Yeung-Jasnow theory of phase ordering 
gives good agreement with simulation results for the correlation function. 

Even though much of the large k behavior of model B can be understood theoretically, a complete theory of phase 
separation, even one that just describes the one-time structure factor, is still needed. 
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In fact, G{x) = vT+~? is the form which results from the YJ theory for model A pO| , pj| . Those authors introduced 
the two-time order-parameter correlation function ||l6|| 



In this case R{t) = [Bt]" with n = 1/2 and B = 4((i— l)/d. Expanding R{ti) and R{t2) in terms of St and r and 
expressing the argument of the arcsin to lowest order in St/t and r/[Bf]"' we get 



ACKNOWLEDGEMENTS 



APPENDIX 




(36) 




(37) 



The expansion for arguments near unity, arcsin(a;) 



n/2 



\/2{l — x), then gives 




(38) 



10 



This correlation function corresponds to Eq. with Ci = 2/7r, C2 = \/ d/8, and G{x) = \/T~fa?. In Ref. |T^] 
the resulting analytic scaling form for the two-time structure factor at large f, given by Eq. (p7|), was obtained 

for model A by a different method utilizing the full form, Eq. (|36|). That approach also yielded a small-f result 
corresponding to the Taylor expansion Eqs. ( |l5|) and the corresponding condition Eq. 
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FIG. 1. Snapshot of a system undergoing spinodal decomposition at the maximum simulation time, r = 4000. All simulations 
reported here were conducted on 512 x 512 lattices. 
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FIG. 2. Central region of the scattering intensity for one realization of the simulation at time r = 4000. The speckling of 
the scattering pattern is apparent, with darker shades indicating brighter speckles. 
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FIG. 3. Time series of the scattering intensity at one wave vector k for one quench to zero temperature. The intensity 

has been normahzod by the timc-dcpcndcnt structure factor averaged over 100 quenches. The dotted hnc is a synthetic 
"Brownian" function constructed to have the same exponential single-time probability density as the scattering intensity and 
an exponential two-time covariance with a characteristic time corresponding to that of the simulation intensity. The persistence 
of the scattering-intensity time series is quite apparent when compared to the synthetic Brownian time series. 



14 



FIG. 4. Schematic illustration of the moving-interface model. The edges of the shaded region represent the interface at 
times ri and T2, and the arrows indicate the directions of motion. The dashed curve is the mean interface. The shaded region 
represents the volume fraction 4>/\, which changes phase between ti and T2. 
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FIG. 5. Probability density of the normalized scattering intensity s{k,T) = I{k,T)/S{k,T), binned for all times and 
0.08 < fe(k) < 0.75. The histogram is exponential for all s observed here, indicating the Gaussian decoupling approxima- 
tion used to equate the squared two-time structure factor and the speckle covariance is a reasonable approximation at the 
one-time level. 
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FIG. 6. Comparison of the two-time structure factor S(k,T\,T2) (squares) and the square-root of the covariance, 
^J'Cav^{k^^\^^^ (circles) for (a) ri — 100, T2 = 200 and (b) ri = 2000, T2 = 4000. These examples serve as a direct 
test of Eq. (^. The agreement is good for small wave vectors, suggesting that two-time correlation functions for real systems 
can be measured via speckle experiments. The arrows indicate the scaling variable z « 4.5 where the Gaussian decoupling for 
two-time correlations breaks down in these simulations. For T2/t\=2, as in this figure, this corresponds to fc'^ri ~ 190. 
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FIG. 7. Contour plot of the scaled speckle two-time correlation function Corr(ti,t2). The correlation decreases as one moves 
away from the line ti =t2, where it is unity by construction. 
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FIG. 8. Characteristic time difference Stc defined by Cov{5tc,t) — (l/2)Cov(0, f ). The hnes are least-squares fits to 
power-law behavior. For t < 50 (dashed hne) the fit gives an exponent 1.02 ± 0.02, which agrees with our expectation for 
small t. For t > 300 (solid line) the least-squares fit gives an exponent 0.66 ± 0.01. This agrees with what is expected from the 
large-t scaling variable z = 5t /t^~" with 1 — n = 2/3, as well as with a recent XIFS experiment [8]. 
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FIG. 9. Scaling of the covariancc vs 5t/t^^^ for several values of f is shown, (a) Linear scale, (b) Linear-log scale, where data 
for the two-time structure factor are also shown (filled symbols). The solid curve is the theoretical form, Eq. (27). Inspection 
of data indicates (27rCi)^ « 32, and numerical fitting to the two-time structure factor data gives nC^ « 0.62. The data collapse 
is good, confirming that there is a single master curve for the speckle covariance in the large-f limit. The deviation of the 
covariance data from the two-time structure factor for z ^ 4.5 indicates the gradual breakdown of the Gaussian decoupling 
approximation for these simulations. 
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FIG. 10. Scaling collapse ol the structure lactor. (a) Scaling of the one-time structure factor for several different times 
presented on a log-log scale vs t = k^r. The slopes of the straight lines correspond to the expected power laws, ^(A:, r) oc k'^ 
at small wave vectors (dashed line) and oc fc"'''"'"^' at large wave vectors (solid line), respectively. The heavy solid curve is 
the theory of Yeung, Oono, and Shinozaki [39]. The inset shows the equivalent scaling collapse of Eq. (3). (b) Scaling of the 

two-time structure factor for fixed ti/t2 — 1/2 and several different ri. This corresponds to a cut through Cov(fc,ri,r2) along 
a line of slope 1/2. The n = 2000 data represented by squares are the same as those represented by squares in Fig. 6(b). 
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FIG. 11. Scaling of the two-time correlation function, (a) For n — 100 and various T2- The zero crossing points remain 
approximately stationary while the amplitude of oscillations decreases monotonically. (b) For several ti/t2 = 1/2. A scaling 
form appears to be approached at later times. The solid curve, the simulation result involving the latest time, is taken as an 
estimate of the scaling function for this ratio of times. 
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FIG. 12. Log- log plot of the autocorrelation vs Tilr\ for several values of n. Least-squares fitting for n = 50, ri > 2000 
gives a slope of 1.49 ± 0.01 which corresponds to A 4.5. This value violates the upper bound A < d conjectured by Fisher 
and Huse [40], but is consistent with the constraint A > (d/2) + 2 derived by Yeung, et al. [41], for conserved systems. 
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